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Abstract 

Orbits of test particles and light rays are an important tool to study the properties of space-time 
metrics. Here we systematically study the properties of the gravitational field of a globally regular 
magnetic monopole in terms of the geodesies of test particles and light. The gravitational field depends 
on two dimensionless parameters, defined as ratios of the characteristic mass scales present. For critical 
values of these parameters the resulting metric coefficients develop a singular behavior, which has profound 
influence on the properties of the resulting space-time and which is clearly reflected in the orbits of the 
test particles and light rays. 

PACS: 

1 Introduction 

To understand the properties of classical solutions of the gravitational field equations it is essential to study 
the orbits of test particles and light rays in these space-times. On the one hand, this is important from an 
observational point of view, since it is only matter and light that is observed and that therefore can give 
insight into a given gravitational field [ ]. The study of the motion of test particles in gravitational fields has 
thus significant practical applications. On the other hand, this study is also important from a fundamental 
point of view, since the motion of matter and light can be used to classify a given space-time, to decode its 
structure and to highlight its characteristics. 

Particles and light have been used since a long time to discuss the properties of solutions of Einstein's 
field equations. All solutions of the geodesic equation in a Schwarzschild gravitational field can be found 
in a seminal paper of Hagihara [2] . With the same mathematical tools one can solve the geodesic equation 
in a Reissner-Nordstrom space-time [3]. The analytic solutions of the geodesic equation in a Kerr and 
Kerr-Newman space-time are also known (see for a survey). Analytic solutions are the starting point for 
approximation methods for the description of real stellar, planetary, comet, asteroid, or satellite trajectories 
(see e.g. [4]). Analytic solutions of the geodesic equation can also serve as test beds for numerical codes 
for the dynamics of binary systems in the extreme stellar mass ratio case and also for the calculation of 
corresponding gravitational wave templates. 

Here we discuss orbits in a static spherically symmetric gravitational field, where the space-time metric 
is expressed in terms of Schwarzschild-like coordinates, 

ds'^ = gttdt^ + grrdr'^ + r'^{dd'^ + sin^ ddip^) . (1) 

We further assume, that the space-time is asymptotically fiat. Well-known examples of such space-times are 
the Schwarzschild solution and the Reissner-Nordstrom solution, both representing black hole solutions with 
an event horizon and a central singularity. In these solutions, the two metric coefficients are not independent, 
but related to each other, gttSrr = —1- At the event horizon g^r — oo, while the curvature invariants remain 
finite. At the origin, in contrast, the curvature invariants diverge, indicating the presence of a classical 
singularity. 

When matter fields other than the Maxwell field are coupled to gravity, not only black hole solutions 
arise but globally regular solutions can appear as well. Interesting examples of such globally regular solutions 
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of the coupled Einstein-matter field equations are boson stars, which form when a scalar field is coupled to 
gravity, or magnetic monopoles, which originate in grand unified theories. These solutions neither possess an 
event horizon nor a central singularity. Their curvature invariants are finite everywhere, and the two metric 
functions no longer satisfy the relation gttgrr — —1, instead gudrr — 

Our interest here focuses on orbits of test particles and light rays in such space-times, and, in particular, in 
space-times emerging in the presence of globally regular magnetic monopoles [5, G] as described in Section 4. 
Magnetic monopoles arise as topological defects in theories which undergo spontaneous symmetry breaking. 
In general, magnetic monopoles exist if the mapping of the vacuum manifold, associated with the symmetry 
breaking, onto the two-sphere is non-trivial. The existence of magnetic monopoles is consequently a generic 
prediction of grand unification. If magnetic monopoles were indeed present in the universe, they would have 
a host of astrophysical and cosmological consequences. One might even think of calculating the form of 
gravitational waves created by scattering processes at such monopole solutions, since these may contribute 
to a stochastic gravitational wave background emerging in the very early universe. 

For magnetic monopoles the Einstein-matter equations depend on two dimensionless parameters, which 
represent ratios of the characteristic mass scales present in the theory. Depending on the values of these 
two parameters, the metric coefficients gu and grr can exhibit an interesting behaviour: the component gu 
then approaches zero within a certain domain, < r < rg, while g^ tends to infinity at the special value 
To, and remains finite elsewhere [(>]. We here demonstrate, that this intriguing behavior is clearly reflected in 
the particle orbits. (When one considers SU(2) Einstein- Yang-Mills theory without Higgs fields, also globally 
regular configurations result [ ]. These solutions are unstable [n], however, and this makes them less attractive 
to study orbits in their vicinity.) 

In Section 2 we recall the general set of equations of motion for test particles and light rays and also 
discuss the effective potential. We then review a number of physical quantities relevant for the interpretation 
of the orbits of the test particles and light rays and discuss the general features of singular limits of the 
metric functions in Section 3. In Section 4 we recall the basic equations and main properties of gravitating 
non-Abelian magnetic monopoles. In Section 5, the main section of the paper, we then present numerical 
results obtained for the possible types of orbits in the gravitational field of a magnetic monopole. We here 
exhibit the trajectories of test particles and light rays together with their proper time resp. affine parameter, 
and we emphasize the peculiar features of the orbits, which arise for almost critical values of the parameters. 
We give our conclusions in Section 6. 

2 Equations of motion and effective potential 
2.1 Equations of motion 

We briefly recall the equations of motion for test particles and light in the general static spherically symmetric 
metric, Eq. (1). Because of spherical symmetry, one can restrict the motion to the equatorial plane. One 
obtains two constants of motion, the specific energy E and the specific angular momentum L (i.e., energy 
and angular momentum per unit mass for a particle). 



E - -gu 



dt 



(2) 




(3) 



where for particles t is the proper time, while for light it represents an affine parameter. Then the only 
dynamical equation left is 



where e = 1 for particles and e = for light. Exploiting energy and angular momentum conservation one 
obtains the equations for </? and t as functions of r 
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Eqs. (4a)-(5b) then give a complete description of the dynamics. From Eq. (4a) it is clear that the motion is 
restricted to the domain E^ + gu (^e + > 0. 

2.2 Effective potential 

When reasoning by analogy to the Newtonian case [ ], a natural definition arises for an effective potential 
of the Schwarzschild space-time from the geodesic equation. For the general asymptotically flat case we are 
discussing here, no such natural definition is available. Therefore we propose a set of conditions which we 
consider to be natural properties of an effective potential, and which in our case are sufficient to give a unique 
definition for an effective potential. 

We thus suggest that an effective potential Ucff for an asymptotically flat space-time should possess the 
following properties: 

1. The effective potential should determine the equation of motion for the r-coordinate: -—^ — 

ar^ dr 

2. The effective potential should obey the boundary condition lim Ucs{r) — 0. 

r — >oo 

The first condition leads to 

^ ^'^''-V-f-C/cff, (6) 



2 \dT J 

where £ is a constant. To fix this constant we decompose the two metric functions gu and grr according to 

^ = -1 - fu{r) , = 1 + frrir) , (7) 

gu{rj grr{r) 

where lim /tt(r) = and lim ,frr{f) = 0, in accordance with asymptotic flatness, and we also exploit the 

r — )-oo r — >oo ' 

second condition. The equation for the radial motion Eq. (4a) then becomes 

^^)' = ^il + fu{r)){l + Urir))E'~^-il + UAr))(^e + 

- + I Uu{r) + frrir) + fu{r)frr{r)) E' ~ y„{r) - (1 + frr{r)) ^ . (8) 

Except for the first term all other terms tend to zero for r — s- cx), so that we can make the unique identification 

£ = , C/cff = -\ ifttir) + frrir) + fu{r)frr{r)) E^ + |/„(r) + (1 + /„(r)) ^ 

2 V gttgrr grr \ J J 

Thus in general the effective potential depends on both constants of motion, on L and on E. 

Since the right hand side of Eq. (8) is quadratic in E, one can recast the equation in the form [10] 

2\dTj 2 2 gtt[r)grr{r) 



and identify the quantity U, 
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where U{r) — > e for r ^ cxd. Although U is sometimes called an effective potential, it is not an effective 
potential in the sense of the above criteria. However, U has the big advantage of being independent of the 
energy E, and from 

E"^ - U{r) = (12) 

one can easily determine the values of r, which mark the turning points and thus give the range of the motion. 
Thus the range of the motion of particles and light is most easily obtained from U and independent of E. 
However, the full features of the geodesies are obtained only from [/cfr and are dependent on E. 



3 Properties of orbits 

For the subsequent discussion of the orbits in the space-time of a magnetic monopole it is instructive to first 
recall a number of physical properties relevant for the description of the motion. These are the proper time, 
the physical distance and the physical velocity. 

We then discuss the limiting features of a metric, when the metric coefficients gu and grr tend to zero or 
diverge. The curvature invariants here indicate the presence of a mere coordinate singularity or of a physical 
singularity. We first address these limits for the Schwarzschild and Reissner-Nordstrom solutions, where 
gttgrr = ^1- Then we turn to the more general case, realized in the space-time of a gravitating magnetic 
monopole, where gttgrr = ~^('')- Here our statements about the limits should be interpreted to only mean 
"very small" and "very large" . The treatment of the limits in the strict mathematical sense will be discussed 
in a subsequent publication [11]. 

These considerations are then used in Section 5 in order to obtain an interpretation of the calculated 
trajectories. They may be also of relevance, when one considers a "splitting" of space-time or a space-time 
"without time" , features which might occur in generalized models of gravity. 

3.1 Proper time 

An important feature of the motion of a particle is the proper time elapsing along its trajectory. Along a 
radial trajectory, e.g., the proper time is given by 




Now if a whole region arises in a space-time, where gtt is very small {gtt — 0), as in the case of the nearly 
critical space-times of a magnetic monopole, then in such a region ^ is very small as well — > 0). This 
means, that a particle does not feel anything like a singularity in such a region, but that it simply needs 
almost no proper time to traverse such a region. 

3.2 Proper distance 

Another relevant physical quantity is the distance as measured from the origin in a regular region of space- 
time, 

I = / Va^d-r ■ (14) 
Jo 

When g^r — > oo at some coordinate value tq, the distance to Tq becomes infinite, if y/grr is not integrable 
at tq. The space-time in the region < r < tq then has infinite extent. 



3.3 Physical velocity 

The velocity dr/dt in Eq. (5b) is the coordinate velocity and not a physical velocity. A physical velocity is 
measured, e.g., by an observer at rest in the given coordinate system. The proper time of such an observer 
at r is given by dT = —^fg^dt, and the physical distance to the origin (in a regular space-time) has been 
defined in Eq. (14). The radial velocity measured by such an observer at rest then is 

dl I grr dr 

(15) 



dT y -gtt dt 
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With the geodesic equation Eq. (5b) this yields for the measured velocity 



The measured velocity is well defined in the allowed regions of the motion, which are restricted by the 
requirement of positivity of the radicand. Clearly, the measured velocity of a particle is always limited by the 
velocity of light. For e = the measured velocity is the light velocity. We note, that the velocity measured 
by an observer at rest does not depend on grr- 



3.4 Curvature invariants 

To differentiate between mere coordinate singularities and physical singularities it is instructive to consider 
the curvature invariants. The simplest of these invariants is the Kretschmann scalar K = R^^p^jR^^^" . For 
the static spherically symmetric metric Eq. (1) the Kretschmann scalar is given by 

K = 




1 9tv 9tt\ 1 9tt 9tt 9rr . 1 9tt\ 9rr\ , 2 / g. 



'^grr 9rr \9tt J 9rr 9tt 9tt 9rr -ig^r \9tt J \9rr J ^^9^X9 



(17) 



In the Schwarzschild and Reissner-Nordstrom solutions, for instance, the Kretschmann scalar remains 
finite at the horizons, whereas it diverges at the classical physical singularities at the origin. 



3.5 Equations of motion for almost singular metric coefficients 

The equations of motion Eqs. (4a)-(5b) depend on the metric coefficients. The equations for dr/dip and dr/dr 
depend on the metric coefficient g^, and on the product gtt9rr, while the equation for dr/dt further contains 
the factor gu- 

To address the effects of almost singular metric coefficients on the orbits of particles and light, we here 
distinguish two cases. In the first case, we consider space-times, where the metric coefficients satisfy the 
relation gtt9rr = ~li such as, e.g., the black hole space-times of the Schwarzschild metric and the Reissner- 
Nordstrom metric. 

In the second case we consider globally regular asympototically flat metrics, where the metric coefficients 
are constrained by the above relation only at spatial infinity, thus gtt9rr — with limr_too^('') — 

1. These metrics include the generic space-time of a magnetic monopole (see section 4). When certain 
critical values of the coupling constants are approached, such a globally regular space-time may, however, 
evolve towards a space-time with singular metric coefficients, possibly associated with physical singularities. 
Studying the orbits of particles and hght is expected to yield a deeper understanding of such limiting space- 
times. 



3.5.1 gttgrr — —1 

When gttgrr — ~1, as, e.g., in the Schwarzschild and Reissner-Nordstrom solutions, the equations of motion 
reduce to 

We now consider a number of effects relevant for the orbits, which arise when grr becomes very small or 
very large, as collected in the following Table: 
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(provided E is large enough) 
* depending on whether Eq. (13) is integrable or not 



When g„ — s- ±oo (g^ — > 0) at some coordinate value tq, a horizon is encountered, as, e.g., in the 
Schwarzschild and Reissner-Nordstrom space-times. The coordinate r then changes from space-like to time- 
like or vice versa, unless the horizon is degenerate as in the case of an extremal Reissner-Nordstrom black 
hole. 

For grr {gu ±oo) on the other hand, a singularity is approached in, e.g., the Schwarzschild and 
Reissner-Nordstrom space-times. When g^,. > 0, motion is allowed only in the region close to the singularity, 
if the energy is large enough to compensate for the large negative second term in the equations (proportional 
to 1/f/rr)- The vanishing of dr (or the finiteness of r) indicates a finite proper time along a geodesic, which 
is a standard criterion for the occurrence of a singularity [ ] . 

3.5.2 gttgrr = 

This case is formulated to be applicable to the motion of particles and light in the gravitational field of a 
magnetic monopole, to be discussed below. The two metric functions are now constrained by the relation 
gttgrr = ^1 only at spatial infinity, while elsewhere their product corresponds to an independent function, 
gttgrr — ~A{r), determined by the Einstein-matter equations. 

We now consider the effects relevant for the orbits, which arise when gu and/or g^r become very small or 
very large: 
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1. When gtt ^ this may 

a) occur at a specific coordinate value vq, and then signal the presence of a horizon, when it is 
associated with g^r ^ ±00. 

b) occur within a whole region of space-time, where > 0. The radial distance then changes 
very fast for small variations of the angle. Therefore a particle crosses such a region on quasi radially 
straight lines. 

2. When gu — > —00, the traversed coordinate range is always finite, while the coordinate velocity tends 
to become large. The latter can be understood since in this case it needs a very small t to generate a 
large proper time. 

3. When grr — > 0, particles need to cover a huge coordinate range in order to move a certain proper 
distance. That means that r has to become very large and, thus, all quantities proportional to dr 
become very large. 

4. When g^r — » 00 at a specific coordinate value rg, 

a) a horizon is encountered, when at the same time gu — >■ 0. The coordinate velocity goes to zero, 
such that the particle appears to be no longer able to move and cross the point vq. 

b) a singularity may arise, leading to a split of space-time. 



4 The gravitating non-Abelian monopole 

In order to be able to discuss the motion of particles and light in the gravitational field of a magnetic 
monopole, we now briefly recall the basic equations and the main features of the monopole space-times. 

4.1 Action 

SU{2) Einstein- Yang-Mills-Higgs theory is described by the action 

^ = / {16^ " (^P.^^^^) - \Tr{D,^D'^^) - ^Tr [{^' - v'f] | V^d'x , (19) 

with curvature scalar R, su{2) field strength tensor 

F^, = df,A, - d,A^ + ie[Af,, A,] , (20) 

gauge potential A^ = A°t°/2, and covariant derivative of the Higgs field $ = <i)'^r° in the adjoint represen- 
tation 

D^^ = d^^ + ie[A^,^] . (21) 

Here G and e denote the gravitational and gauge coupling constants, respectively, v denotes the vacuum 
expectation value of the Higgs field, and A represents the strength of the Higgs self-coupling. 

The nonzero vacuum expectation value of the Higgs field breaks the non-Abelian SU (2) gauge symmetry 
to an Abelian U{1) symmetry. The particle spectrum of the theory then consists of a massless photon, two 
massive vector bosons of mass Mw = ev, and a Higgs field of mass Mh = V^Xv. In the limit A — > the 
Higgs potential vanishes, and the Higgs field becomes massless. 

4.2 General equations of motion 

Variation of the action Eq. (19) with respect to the metric g'^'^ leads to the Einstein equations 



G^,y — R^u — -^^tivR — SnGTf^i, (22) 



with stress-energy tensor 



T^, = 2Tr(F^,F,^5"'5- Js^.F^^F"^) 

+ Tri^-D^^D,^ - ^g^.D^^D'^^) - ^g^,Tr{^^ - . (23) 
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Variation with respect to the gauge potential and the Higgs field <& leads to the matter field equations, 

^ =Z?„(^F^^) - \ie[<^>, i?"*] = , (24) 



-^Df,{y/^D^'<^>) + X{^^ ~v^)<i> ^0 . (25) 

4.3 Ansatze 

To obtain a static spherically symmetric magnetic monopole with unit magnetic charge we parametrize the 
metric in Schwarzschild-like coordinates [5, (i] 

ds^ = gtt{r)df + grr{r)dr^ + r^{de^ + sin^ Odif)'^) , (26) 

and introduce the mass function m(r) 

grAr)=(l-'^y . (27) 
For the gauge potential and the Higgs field we employ the spherically symmetric Ansatze ['i, (>] 

A^dx'' = ^—^^ (r^dO - Te sin Odip) , (28) 

and 

$ = vH{r)Tr , (29) 

where the su{2) matrices r^, Tg and are defined as products of the spherical spatial unit vectors with the 
vector of Pauli matrices t° . 

4.4 Dimensionless quantities 

We now introduce the dimensionless coupling constants a and /3 as ratios of the mass scales present in the 
theory, 

, ^AttM^ ^ ^ 1 _ A 

" -7^1;-^^^'' ' ^-2M^-^' 

and Mpi = 1/\/G is the Planck mass. We further pass to dimensionless coordinates, evr r, and a 
dimensionless mass function, evm m. 

4.5 Monopole equations of motion 

The tt and rr components of the Einstein equations then yield equations for the mass function m and for the 
product gttgrr, 



m' ^ a'i K'' + -r^H'' - —K'^ - mvH'^ + + h^k^ + (^r\H^ - 1)M , (31) 



and 



{gttgrr)' = 2a\\^^ + H'''j{gugrr) , (32) 
where the prime indicates the derivative with respect to r. For the matter functions we obtain the equations 

yf-gugrr \ grr J \ J 
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and 



^ (" ^V gugrrH' y ^ ^ ^2i^2 _^ p\\H^ _ 1)) . (34) 

This set of equations depends only on the dimensionless couphng constants a and (3. Since no analytic 
solutions to these field equations are known, we use numerical solutions for the calculation of the test particle 
orbits. 

4.6 Boundary conditions 

To obtain globally regular particle-like solutions we require at the origin the boundary conditions 

m(0) = , K{Q) = 1 , iJ(0) = . (35) 

Asymptotic flatness implies that the metric functions gu and g^r approach constants at infinity. We adopt 

gtt{oo)grr{oo) = -1 . (36) 

m(oo) represents the dimensionless mass of the monopole solutions. 
The matter functions satisfy asymptotically 

K{oo) = , H{oo) = 1 . (37) 



4.7 Embedded Reissner-Nordstrom solutions 

A special black hole solution of the set of coupled equations is the embedded Reissner-Nordstrom solution 
with unit magnetic charge, 

m(r) = TOoo - — , gtt{r)grr{r) ^ -I , (38) 
2r 

K{r) = , H{r) = 1 . (39) 
The extremal Reissner-Nordstrom solution with event horizon radius rjj satisfies 

= moo = a . (40) 



4.8 Properties of magnetic monopole solutions 

The coupling to gravity has a significant effect on the magnetic monopole solutions present in flat space . 
When the coupling constant a is increased from zero, a branch of gravitating monopole solutions emerges 
smoothly from the flat space 't Hooft-Polyakov monopole solution. This branch of gravitating monopole 
solutions extends up to a maximal value amaxi beyond which gravity becomes too strong for regular monopole 
solutions to persist [fi]. 

For vanishing coupling constant /?, this first gravitating monopole branch merges with a second branch 
at Q!max, which extends slightly backwards, up to a critical value acr of the coupling constant. At acr this 
second branch of gravitating monopole solutions bifurcates with the branch of extremal Reissner-Nordstrom 
solutions. 

In particular, at rp = a double zero of the metric function \/grr appears. But does not correspond 
to a degenerate horizon. The Kretschmann scalar diverges there. The exterior space-time of the solution, 
however, corresponds to the one of an extremal Reissner-Nordstrom black hole with unit magnetic charge 
[()]. The interior space-time retains regularity at the center, due to the influence of the non-Abelian flelds 
present. 

As (3 increases, the second branch decreases in size, until at a certain value of [3 the maximal value ctmax 
and the critical value acr coincide [(>]. For a considerably larger value of (3 another interesting phenomenon 
arises. The metric function 1/grr develops a second minimum, and the double zero now arises at a value 
r* < ro, where r* does correspond to a degenerate horizon. The critical solution then corresponds to an 
extremal black hole with non-Abelian hair and with a mass less than that of the corresponding extremal 
Reissner-Nordstrom solution [13]. 
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5 Motion in the space— time of a gravitating monopole 



This section constitutes the central part of the work. Here we present and discuss the possible types of 
orbits in the space-time of a gravitating monopole. The monopole solutions depend on the dimensionless 
parameters a and p. The parameters then enter via the respective metric functions in the equations of motion 
and determine via the potentials the allowed trajectories of particles and light. 

We first address the potential U and the effective potential UeS in these monopole space-times, considering, 
in particular, their expansion close to the origin. We then discuss the orbits of particles and light obtained 
numerically in these space-times. Here we exhibit generic examples of the possible types of motion, and we 
discuss motion in space-times very close to critical monopole solutions. 

One of the interesting features of the motion in the space-time of a monopole is the capture of light rays 
by the source when a is big: here light is found to move on bound geodesies. Large a space-times also admit 
two bound regions for particles to move in. 



5.1 Potentials 

The space-time of a generic non-Abclian magnetic monopole is globally regular and asymptotically flat. 
Only when critical values of the parameters are approached the space-time may evolve coordinate or physical 
singularities. 

Asymptotic flatness implies appropriate boundary conditions for the metric functions, which are reflected 
in the fall-off of the functions f^r and fu, introduced in Eq. (7) to define the effective potential. From Eq. (27) 
we identify f„ for the monopole space-time, 

/..W--"^ (41) 
r 

which vanishes for r ^ oo. From the monopole boundary condition Eq. (36) we also infer that fu — >■ for 
r — > oo. Thus the effective potential C/eff has the required asymptotic behaviour. 

The generic monopole space-time is globally regular, i.e., there is, in particular, no singularity present at 
the origin, but the space-time is smooth there. This is seen in the expansions of the metric functions at the 
origin, since for small r one finds [6] 

/„(r) = -caV2 + 0(r4), fu{r) = c' + 0{r) , (42) 

where c and c' are constants. It is also reflected in the smooth behaviour of the curvature invariants close to 
the origin. 

Close to the origin, the potentials then have the following expansion 

C/cff = ^+d+0{r), (43) 
2r^ 

U = f^ + ^+d" + 0{r), (44) 
2r^ r 

where d, d' , and d" are constants. Therefore, for r — *■ a repulsive angular momentum barrier is always 
present in the effective potential, unless the motion is purely radial. In orbits with finite angular momentum 
particles or light rays can never reach the origin, r = 0. 

As in the Newtonian case, the angular momentum barrier here dominates close to origin. But while the 
angular momentum barrier increases with increasing L, its influence is modulated by the metric function gu, 
so that in regions with very small \gtt\ the influence of the angular momentum barrier is reduced. This is 
seen in Fig. 1, where the potential U and the effective potential Ucs are shown for {a, /3) = (0.85,8.282558). 



5.2 Orbits for /3 = 

We now discuss the possible orbits for particles and light rays in monopole space-times at vanishing (3. We 
begin with a set of typical solutions, obtained for (a,/3) = (0.25,0). As a is increased a maximal value Ofmax 
is reached, which still allows for the existence of a globally regular monopole solution. Beyond amax only 
black hole solutions exist. We exhibit a set of typical orbits present at amax = 1.403. 
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Figure 1: The potential U (a) and the effective potential UcS at energy E 
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parameters (o,/3) = (0.85,8.282558) and several angular momenta L (see subsection 5.3.2). 



(b) versus the radial coordinate r for the 



When (3 — (or small), two regular monopole solutions exist in the range aa < a < amax- For P — 0, 
the critical value of a is determined to be (within numerical accuracy) — 1.385853. It is called critical, 
because as a — > ctcr the minimum of the metric function 1/grr decreases, and reaches zero at a certain value 
To of the radial coordinate in the limit. The region of space-time with r > rg of the critical solution then 
corresponds to the exterior space-time of an extremal Reissner-Nordstrom black hole. We exhibit a set of 
orbits at ttcr = 1.385853, i.e., for an almost critical space-time. 

In the following we always present figures for the metric functions —gtt and l/gm the curvature invariant 
K, and the potential U, and exhibit various typical and special orbits obtained for given values of the angular 
momentum L and the energy E. 



5.2.1 Orbits at a = 0.25 

We begin with the discussion of the generic case of a monopole space-time, as exhibited in Fig. 2 for (a, (3) = 
(0.25, 0). The metric coefficient —gtt has a finite value at the origin and rises monotonically to its asymptotic 
value, — gtt (oo) = 1. The metric coefficient l/grr assumes the value 1/grr = 1 both at the origin and 
asymptotically, and exhibits a minimum at some value of the radial coordinate, tq. The Krctschmann scalar 
is finite everywhere and small. 

The potential U shows, that the orbits in such a space-time have the general structure of typical Newtonian 
orbits: there is one circular orbit corresponding to the minimum of U, and there are bound orbits as well as 
scattering states. 

To exhibit the possible types of orbits, we choose as a representative value for the angular momentum, 
L = 0.5. As seen in Fig. 3, the orbits indeed show the qualitative structure inferred from U. The main 
difference to the typical Newtonian orbits is the occurrence of a perihelion (perimonopolion) shift for bound 
orbits, showing that the underlying potential differs from the Newtonian potential of a point mass. The 
elapsed proper time along the orbits is also indicated. 



5.2.2 Orbits at a^ax = 1.403 

We now consider the orbits for the maximal value of a, i.e., for the parameters (amaxi/3) = (1.403,0). 
Again, the metric coefficient —gtt is a monotonically rising function, and the metric coefficient 1/grr exhibits 
a minimum at some value of the radial coordinate, tq, although the minimum is much deeper now. The 
Kretschmann scalar is still finite everywhere and not too large. 

The most interesting feature here is that the potential U can now possess two minima and a maximum. 
In such a case, bound orbits are present in two distinct regions of space. We exhibit typical examples of 
such bound orbits in Fig. 5 for an angular momentum of i = 4.3. The first three figures (a), (b) and (c) 
show orbits in the inner region (associated with the inner minimum), for a low energy value, E = -^0.2, for 
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Figure 2; The metric functions —gtt (a) and l/g^r (b), the potential U (c) and the curvature invariant K (d) versus the radial 
coordinate r for (a,/3) = (0.25,0). 




Figure 3: Particle orbits r{if) for {a, fS) = (0.25,0). The dots indicate units of elapsed proper time. 



a higher energy vahic, E = \/(L5, and for the (all but) highest value possible for a bound orbit in the inner 
region, E = VO-908. 

For very low energies the orbits of course approach a circular orbit. As the energy increases, the orbits 
have a deformed shape, and show a very pronounced perihelion shift ((a) and (b)). When the limiting energy 
value E'max for bound orbital motion in the inner region is approached, the associated unstable circular orbit 
becomes apparent, since it forms the outer envelope of the trajectory (c). 

In the energy interval between the outer minimum and the maximum, there are orbits in the inner region 
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Figure 4: The metric functions —gtt (a) and 1/grr (b), the potential U (c) and the curvature invariant K (d) versus the radial 
coordinate r for (amaxi/3) = (1-403,0). 



(c) as well as in the outer region (d) . Here the unstable circular orbit associated with the maximum becomes 
apparent as the inner boundary for the orbital motion. For still higher energies, in the range i?max < i? < 1 
one observes quasi-elliptic orbits (e). In the outer region these look like ellipses. When they get close to the 
inner region, the unstable circular orbit is again reflected in the motion, but it no longer forms a boundary. 
Instead the motion proceeds further inwards, orbits close to the center and then proceeds outwards again, as 
highlighted in the inset of (e). 

For energies i? > 1, we observe scattering states (f). The particle then encircles the central region several 
times (twice in (f)) before escaping into infinity again. 

5.2.3 Orbits at = 1.38585 

We now demonstrate, how the monopole space-time evolves as the critical solution is approached. As a — > acr 
the minimum of the metric function 1/grr decreases and reaches zero in the limit. This decrease towards zero 
is seen in Fig. 6. Note, that acr is here determined only within a certain numerical accuracy, and therefore 
the final solution obtained is only almost critical: 1 / grr does not yet fully reach zero at some coordinate value 
ro, but it is extremely small there. 

The true limiting solution then consists of two parts: the interior region < r < rg and the exterior 
region < r < oo. Since the limiting exterior solution corresponds to an extremal Reissner-Nordstrom 
space-time, rg = a must hold for the limiting solution, according to Eq. (40). The limiting interior solution 
in a non-Abelian solution and regular at the origin. 

As a — > ttcr the metric function —gu becomes increasingly small in the inner region, Q < r < r^, tending 
to zero in the limit. Close to rg, however, —gtt rises very steeply to assume its asymptotic value of one, as 
imposed by the boundary conditions. In the limit a — *■ acr, —gtt appears to become singular at rg. 

The emergence of a singularity at ro is also indicated by the Kretschmann scalar of the monopole space- 
time, which appears to diverge in the limit. In contrast, the Kretschmann scalar of the associated extremal 
Reissner-Nordstrom black hole remains finite at its degenerate horizon tq, but diverges for r ^ 0, as seen in 
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Figure 5: Particle orbits r(ip) for («max,/9) = (1.403,0). The crosses indicate units of elapsed proper time. 



Fig. 6. 

Let us now consider the motion of particles in the almost critical space-time, with parameters (a, P) — 
(1.38585,0), and exemplify the intriguing effects of this almost critical space-time on the orbital motion. 
We study the orbits for a representative value of the angular momentum, L = 4.3. As shown in Fig. 6, the 
potential U then again possesses two minima. We therefore expect the same general types of orbits as in the 
previous case. 

We display a set of characteristic orbits in Fig. 7. The bound orbits in the inner region (a) and (b) 
(connected to the inner minimum of the potential) are now hugely distorted as compared to ordinary bound 
orbits. In the energy interval between the outer minimum and the maximum, bound orbits are present in the 
inner (c) as well as in the outer region (d). Again the unstable circular orbit associated with the maximum 
becomes apparent as the outer envelope of the allowed motion in the inner region (c) and as the inner 
boundary for the trajectories in the outer region (d). For still higher energies, in the range -Emax < -E < 1 
we again observe quasi-elliptic orbits with additional inner loops circling the center (e), while for energies 
-B > 1, we again observe scattering states (f). 

The most surprising feature is the deformation of the orbits in the inner region. As seen in Fig. 6, the 
orbits evolve smoothly until they reach the vicinity of the critical value of the radial coordinate, rg. There 
they abruptly change direction and approach in almost straight radial lines the center. Very close to the 
center they are reflected (passing partial arcs), and then move outwards again in almost straight radial lines, 
until they reach the vicinity of rg, where they again abruptly change direction. 

This intriguing pattern of movement can be understood as follows: for r < rg the effective potential very 
quickly tends to extremely large negative values, as seen in Fig. 6. Therefore in this region a particle is 
vigorously attracted towards the origin. This huge attraction leads to orbits which are almost straight radial 
lines. (Note, that this interpretation is related to the first column of the table in subsection 3.5.2.) The 
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Figure 6: The metric functions —gtt (a) and 1/grr (t>), the curvature invariant K (c) and (d), the potential U (e), and the effective 
potential U^ff (f), versus the radial coordinate r for (ocr, /3) = (1.38585, 0). (RN denotes the extremal Reissner-Nordstrom black 
hole.) 



deflection close to the center is caused by the potential barrier there. 

Interesting is also a glance at the proper time of a particle along these almost straight radial parts of the 
orbits. Since —gtt is extremely small, a particle needs almost no proper time to traverse the inner region, as 
exemphfied in Fig. 6 (and also discussed in subsection 3.5.2). 

Let us finally consider these orbits in terms of the proper distance / instead of the radial coordinate r. In 
Fig. 8 we exhibit the dependence l{r) for a sequence of solutions approaching the critical solution. The figure 
demonstrates, that the space-time develops a throat, as a ^ acr- In the limit, the throat becomes infinitely 
long [()]. For the almost critical value a = 1.385853, I grows already by a factor of 20 in the vicinity of Tq. 

The rapid growth of I in the vicinity of ro has significant infiuence on the shape of the orbits. The abrupt 
directional changes seen in the orbits r{ip) in the vicinity of rp are smoothed out, when l{ip) is considered 
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Figure 7: Particle orbits r{ip) for («cr,/3) = (1.385853,0). The dots indicate units of elapsed proper time. 



instead, since now the orbits are seen to move within the long throat of the space-time. This is demonstrated 
in Fig. 9 for the inner orbits of Fig. 6. 



5.3 Orbits for /3 > 

Let us now consider the effect of a finite vakie of the parameter (3 and thus the Higgs mass on the monopole 
space-times and consequently on the orbits in these space-times. 

As f3 increases amax decreases. At the same time, and ctmax approach each other, until they merge. 
The critical space-time then arises at the maximal possible value of a for the given value of (3. We exhibit 
such an almost critical space-time at {a, (3) = (0.9,4.49088) and study its orbits. 

As (3 increases further another interesting phenomenon arises: l/g^r develops two minima. We study the 
evolution of such space-times towards the corresponding critical space-time at {a, (3) — (0.85,8.282558). We 
discuss the effects on the orbits of particles, and we exhibit orbits of light rays. 



5.3.1 Orbits at o; = 0.9, /3 = 4.49088 

We now address motion in an almost critical space-time at a finite, but still small value of P, namely at 
P — 4.49088. This value of /3 is obtained, by fixing a = 0.9 and then increasing /3, until a = 0.9 becomes the 
maximal and at the same time critical value, oimax = Ocr = 0.9, i.e., for this specific value f3 = 4.49088 no 
regular solutions exist beyond a = 0.9. 

We exhibit the metric coefficients of this almost critical space-time in Fig. 10. The evolution of the 
monopole space-time towards the critical solution is similar the one discussed in the previous case, (a, f3) = 
(1.38585,0). As a the minimum of the metric function 1/grr decreases and reaches zero in the limit. 
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Figure 8: The measured spatial distance / versus the radial coordinate r. As a approaches a^rit a- throat develops at r = tq. 




(a) (b) (c) 



Figure 9: Particle orbits l{'p) for (Q!cr,/3) = (1.385853,0). The dots indicate units of elapsed proper time. 

The true limiting solution then consists of the non-Abelian interior solution in the region < r < tq, and the 
extremal Reissner-Nordstrom solution with degenerate horizon at rp = a in the exterior region tq < r < oo. 

Again, the metric function —gtt becomes increasingly small in the inner region, 0<r<ro, asa— > acr, 
and rises very steeply in the vicinity of tq. The main difference to the previous case is, that —gtt is no 
longer monotonic, but that it has a minimum close to tq, before the steep rise. (Note, that this minimum 
foreshadows the appearance of a second minimum in grr for still larger values of /3, discussed below.) In the 
limit a aa, ~gtt again appears to become singular at tq, as reflected in Fig. 10 by the Kretschmann scalar. 

Considering the motion of particles in the almost critical space-time, we again select for the angular 
momentum the value L — 4.3. The potential U then again possesses two minima, giving rise to similar orbits 
as observed in the previous case. 

We display a set of characteristic orbits in Fig. 11. At very low energy (a) we observe smooth bound 
orbits in the inner region, which are limited by a circle at tq, that appears to form an outer envelope for the 
motion. As the energy increases, the bound orbits in the inner region (b) and (c), are again hugely distorted, 
exhibiting smooth segments alternating with almost straight radial lines. 

The energy for the bound orbit in (b) is chosen, so that the motion is still limited by the circle at tq, 
but that the motion now either (almost) proceeds on segments of that circle or that it proceeds on almost 
straight radial lines to and from the center (with some reflection there). 

The particle motion is again understood by examining the effective potential. As before, the abrupt 
changes in the motion followed by the straight radial line segments are caused by the extreme decrease of the 
effective potential in the region r < tq, leading to vigorous attraction for the particle towards the origin. At 
the same time the particle needs almost no proper time to traverse the inner region. 
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Figure 10: The metric functions —gtt (a) and 1/grr (b), the potential U (c), and the curvature invariant K (d) versus the radial 
coordinate r for (a, /3) = (0.9,4.49088). (RN denotes the extremal Reissner— Nordstrom black hole.) 



Fig. 10 also exhibits a bound orbit in the outer region (d), and a scattering state (e). This scattering orbit 
reflects the maximum of the potential, when the particle traverses the inner loops circling the center, as well 
as the critical value rp, since the particle passes also the associated straight radial line segments close to the 
center. 



5.3.2 Orbits at a = 0.85, /3 = 8.282558 

As our last example for motion in the space-time of a gravitating monopole, we discuss a case with interme- 
diate Higgs mass, choosing the parameter set = (0.85,8.282558). Here a new phenomenon arises: the 
limiting space-time no longer corresponds to an extremal Reissner-Nordstrom solution in the exterior, but it 
retains non-Abelian fields there, which affect the features of the space-time [ ]. 

Let us first consider, how the monopole space-time evolves in this case as the critical solution is ap- 
proached. As in the last example, we fix a and increase /?, until the almost critical space-time is reached. 
As /3 — > Per (or analogously as a ^ act) the minimum of the metric function l/g^r close to — a decreases, 
but it does not reach zero. Instead, at a certain value of the parameter a second minimum appears at a value 
r* < rg. As the parameter is further evolved towards its critical value, it is this new minimum at r* , which 
reaches zero in the limit, as anticipated from Fig. 12. The old minimum at vq is no longer affected by the 
final approach towards the critical solution. 

The true limiting solution thus divides space-time in a different way and consists of the interior region 
< r < r* and the exterior region r* < r < oo, where the exterior region may be subdivided into the 
intermediate region r* < r < rg and the outer region rg < r < oo. Clearly, the metric of the limiting solution 
differs now from the metric of an extremal Reissner-Nordstrom black hole also in the exterior region of the 
limiting space-time r* < r < oo. Though for rg < r < cx3 it is still close to an extremal Reissner-Nordstrom 
space-time for the chosen parameters. 
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Figure 11: Particle orbits (a)-{e) r{ip) and (f) l{ip) for (a,/3) = (0.9,4.49088). The dots indicate units of elapsed proper time. 



As in the case above with a single minimum at tq, the metric function —gu becomes very small in the 
region < r < tq, when (3 Per- However, unlike that case, it decreases distinctly faster at r* than in the 
overall region < r < tq. Thus —gtt develops a minimum at r* quite analogous to the minimum of 1/grr- 
As a consequence, the product of —gu and grr changes only slightly in the vicinity of r* , while, in contrast, 
~gttgrr changes very rapidly in the vicinity of j'q. 

The implication for the curvature of the space-time is, that there appears no singularity at r* . The 
Kretschmann scalar remains perfectly smooth there, as seen in Fig. 12. On the other hand, the Kretschmann 
scalar still becomes very large at rg, but it will not diverge there in the limit. The critical space-time is thus 
expected to correspond to a space-time with a degenerate black hole horizon at r*, and with no singularities 
at either or at the origin. 

Let us now address the orbits of particles and also of light rays in this almost limiting space-time. We 
again choose for the angular momentum the previous value L = 4.3. The corresponding potential U is seen 
in Fig. 12. 

Time-like geodesies 

We display a set of characteristic particle orbits in Fig. 13. Since r* is a minimum of gu the potential U 
also exhibits a minimum at r* (the term is negligible here). For a very small minimal energy this implies 
an almost circular motion at r* , exhibited in (a). The associated effective potential Ucs (Fig- 12 (f)) also 
assumes its minimal value at r* . 

As the energy is increased, the effective potential changes drastically and rapidly assumes huge negative 
values in the region bounded by vq (and not only the region bounded by r*). Thus tq still plays a most 
significant role for the orbits. For instance, the energy for the bound orbit in (b) is again chosen, such that 
the motion is limited by the circle at rg. The motion then proceeds either on segments of that circle or on 
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Figure 12: The metric functions —gtt (a) and l/flrr (b), the potential U (c), and the curvature invaxiajit K (d), the effective 
potential C/eff (c) and (f), versus the radial coordinate r for {a, 13) = (0.85,8.282558). (RN denotes the extremal Reissner- 
Nordstrom black hole.) 
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Figure 13: Particle orbits r{(p) for {a, fS) = (0.85,8.282558). The dots indicate units of elapsed proper time. 
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Figure 14: The potential U for null geodesies versus the radial coordinate r for (a,/3) = (0.85,8.282558), as well as for 
(a, 13) = (1.403,0) (subsection 5.2.2). 
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Figure 15: Null geodesies r((^) for (a,/3) = (0.85,8.282558). 



almost straight radial lines to and from the center (with some reflection there). 

At somewhat higher energies the bound orbits in the inner region (c) and (e) again exhibit smooth segments 
which alternate with almost straight radial lines to and from the center, where the abrupt directional changes 
occur at tq (and at the center). This remarkable role of ro is preserved also in the weakly bound orbits (g) 
and in the escape orbits (h) and (i). As discussed above, a particle needs very little proper time to traverse 
the inner region. 

The bound orbits in the outer region (d) and (f), of course, do not penetrate far enough inside to experience 
any influence of rg. These trajectories are ordinary quasi-elliptic orbits with a perihelion shift. 
Null geodesies 

We have also considered the null geodesies in this space-time. The potential U for light rays (e = 0) is 
shown in Fig. 14, while a set of orbits is exhibited in Fig. 15. The null geodesies are interesting due to the 
fact, that light rays can follow bound geodesies, as evident from the form of the potential U. Thus light can 
be captured by the source, but does not fall into any singularity, since the space-time is regular. 

The shape of the null geodesies is similar to the shape of the time-like geodesies, as seen in Fig. 15. The 
bound orbits (a) again exhibit smooth segments which alternate with almost straight radial lines to and from 
the center, with abrupt directional changes at ro (and at the center), related to the fact that the effective 
potential becomes very large and negative in the region within the radius rg. This peculiar effect on the 
motion at and within radius rg is also seen in the scattering orbits at high enough energy (c) . The scattering 
orbits at lower energy (b) are deflected before they can experience any influence of rg. 
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6 Conclusions and outlook 



An accurate analysis and true interpretation of a given gravitational field can only be obtained through the 
exploration of the geodesies of particles and light rays in this space-time. Such an investigation is particularly 
important in regions, where the space-time metric exhibits unusual behavior, as for example in the vicinity 
of a black hole event horizon. 

In this paper we have discussed the motion of particles and light rays in the gravitational field of a 
magnetic monopole which is characterized by two dimensionless constants, a and (3, where a signifies the 
relative strength of gravity, while /3 is the ratio of the Higgs boson to the vector boson mass in the non-Abelian 
gauge theory. 

The space-time of a non-Abelian monopole is globally regular. Its curvature invariants are finite ev- 
erywhere. However, globally regular monopoles cannot exist for large values of the gravitational coupling 
strength. Simple dimensional reasoning shows, that black holes should form, when a reaches values on the 
order of one [(>]. Indeed, at critical values of the parameters (acr,/3cr) the space-time changes dramatically 
and the metric becomes singular. 

As the critical space-time is approached for small values of the metric coefficient 1/grr tends to zero 
at tq = acr, and the metric coefficient gtt tends to zero in the interval < r < tq. The effects on the particle 
orbits are then astounding. The orbits traversing the region, where tq is located, are hugely distorted as 
compared to ordinary bound orbits. While the orbits evolve smoothly until they reach the vicinity of rg, 
they then change abruptly direction and approach the center in almost straight radial lines. There they are 
reflected, and move in almost straight radial lines out to vq, where they change abruptly direction again and 
evolve smoothly further. This intriguing pattern of movement is caused by the steep drop of the effective 
potential for r < rg, leading to a vigorous attraction of particles or light rays in this region. Also, particles 
need almost no proper time to traverse this region (r < rg). 

For larger values of /3, the critical space-time is distinctly different [1-^]. The metric coef&cient l/grr 
acquires a second minimum at r* < rg, which tends to zero in the limit, while the minimum at rg remains 
finite. Still it is rg, and the steep drop of the effective potential associated with it, which strongly dominates 
the orbits in the interior, causing abrupt directional changes and vigorous attraction towards the center, along 
with almost no lapse of proper time. 

We note, that this kind of research may be extended to metrics associated with stationary axially symmet- 
ric solutions of the Einstein field equations. There are interesting examples where counterrotating horizons 
appear [14] or negative horizon masses [lO]. Such unusual features of various space-times may be explored 
best through the study of the motion of particles and light. This may be extended to the study of particles 
with spin [; ;;]. 

Another extension of the present considerations is the study of the motion of satellites and stars in the 
gravitational field of a modified gravitational theory. This may have applications to the Pioneer anomaly, 
the fiyby anomaly, or the increase of the astronomical unit [ ] , which are all problems which are unresolved 
within standard general relativity. 
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